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TECHNICAL NOTE 5658 


ON PANEL FLUTTER AND DIVERGENCE OF INFINITELY 
LONG UNSTIFFENED AND RING-STIFFENED 
THIN-WALLED CIRCULAR CYLINDERS 
By Robert W. Leonard and John M. Hedgepeth 

SUMMARY 


A preliminary theoretical investigation of the panel flutter and 
divergence of infinitely long, unstiffened and ring-stiffened thin-walled 
circular cylinders is described. Linearized unsteady potential -flow 
theory is utilized in conjunction with Donnell’s cylinder theory to 
obtain equilibrium equations for panel flutter. Where necessary, a sim- 
plified version of Fliigge's cylinder theory is used to obtain greater 
accuracy. By applying Nyquist diagram techniques, analytical criteria 
for the location of stability boundaries are derived. A limited number 
of computed results are presented. 


INTRODUCTIQN 


Although considerable effort has been expended in studying the flut- 
ter and divergence of thin, flat panels exposed to an airstream (see, for 
example, refs. 1 to 9) , little is known of the importance of similar aero- 
elastic phenomena in the design of thin-walled cylindrical missile bodies 
or of other aircraft components where thin, curved panels are used. The 
purpose of this report is to describe a preliminary theoretical investiga- 
tion of the aeroelastic stability of such configurations. Analytical 
criteria for the determination of panel flutter and panel-divergence 
boundaries for infinitely long, unstiffened and ring-stiffened thin-walled 
circular cylinders are presented along with a limited number of computed 
results . 






HACA W 5658 


SYMBOIfi 


distance between ring stiffeners 

amplitude of mth term in expansion for lateral motion of 
ring-stiffened cylinder 

dummy variable 

speed of sound in air 

speed of sound in fluid inside cylinder 


speed of sound in cylinder material, 



plate flexural stiffness per unit length. 


Et' 


12(1 - |i 2 ) 


Young's modulus 

base of natural system of logarithms 

outside and inside air-force functions defined by equations (5) 
and (6), respectively 

outside- and inside air-force functions for vibrating ring (see 
eqs. (3*0 and (35 ) > respectively) 


Hankel functions of first and second kind, respectively, 
of order n 

thickness of cylinder wall 
imaginary part 

modified Bessel function of first kind of order n 


Bessel function of first' kind of order n 


integer 

modified Bessel function of second kind of order n 


dimensionless frequency of harmonic vibration, for unstiff 

ened cylinder and — for ring-stiffened cylinder 
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L 

L 

l 

l 


M 

Mi 

m 


N 


n 

p 0^ p l 


P 0^ P 1 


R 

Re 

r 

t 

w 

w 


X 

X I 


resonant frequencies of internal air-force function for 
unstiffened cylinder (see eq. (22)) 

function defined by equation (8) 

function defined by equation ( 9 ) 

outward lift force on cylinder vail 

amplitude of harmonically varying outward lift force for 
unstiffened cylinder 

amplitude of muh term in expansion for outward lift force for 
stiffened cylinder 

Mach number of flow along vibrating cylinder 

Mach number of external or internal flow along stationary 
cylinder 

integer, number of longitudinal half-waves in each bay of 
ring-stiffened cylinder 

function defined by equation (44) 

integer, number of full waves around circumference of cylinder 
forces exerted on cylinder wall by ring stiffeners 

amplitudes of harmonically varying reaction forces exerted 
by ring stiffeners 

radius of cylinder 

real part 

radial coordinate 

time 

lateral deflection of cylinder wall, positive outward 

amplitude of lateral motion of cylinder wall for unstiffened 
cylinder 

longitudinal coordinate for vibrating cylinder 
longitudinal coordinate for stationary cylinder 
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z 

argument of Bessel functions for supersonic 

relative flow, 



f v t v/l 2 - 1 for unstiffened cylinder and 
for stiffened cylinder 

Met J \ 2 - 1 

i 

a 

aspect-ratio parameter, Rn/a 


y' 

Pm 

function defined after equation (39) 



6(x) 

Dirac delta function, S(x) = 0 for x / Oj 

Poo 

1 5(x)dx = 1 



(J-oo 


e damping coefficient 

£ argument of Bessel functions for subsonic relative flow, 

j v I \j 1 - |2 for unstiffened cylinder and |mja \j 1 - I** 
for stiffened cylinder 

0 angular circumferential coordinate 

p Poisson's ratio 

v wave number of longitudinal waves in unstiffened cylinder 

v dimensionless wave number, Rv 

£ dummy variable 

p mass density of air 

mass density of fluid inside cylinder 

p_ mass density of cylinder material 

a Q) a x midplane stresses in circumferential and longitudinal 
c s ‘ directions, positive in tension 

cr 0 ,a x dimensionless midplane stresses, ctq/E and cr x /E, respectively 

ft parameter defined after equation (8) 

fijjj parameter defined after equation (39) 

a) frequency of harmonic vibration 
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(2L + i JSLf 

VSx 2 - R 2 S0 2 / 
operator 


cr 

critical value 

min 

minimum value 

max 

maximum value 


Primes are used to indicate differentiation with respect to complete 
argument. Subscript notation is used to denote partial differentiation. 


METHOD OF APPROACH 


The configuration under consideration consists of a thin-walled 
unstiffened or ring-stiffened circular cylinder extending to infinity in 
'■* the positive and negative x-directions . (See fig. 1.) The cylinder is 

filled with a stationary fluid and is surrounded by air flowing in the 
positive x-direction at a Mach number M. The effects of midplane ten- 
sile stresses in both the circumferential and longitudinal directions 
and of a small amount of structural damping are taken into account. 

For simplicity in the analysis, it is assumed that the deformations 
of the cylinder walls can in most cases be adequately described by 
Donnell's differential equation. (See ref. 10.) It is kept in mind, 
however, that the validity of Donnell’s theory is limited to cases in 
which there are a large number of circumferential waves; where this con- 
dition is violated, a simplified version of Fliigge's cylinder theory 
(see refs. 11 and 12) is employed to improve the accuracy of the results. 

The problem to be considered is the determination of those combina- 
tions of the parameters characterizing the cylinder and its environment 
that correspond to the boundary between states of stable and unstable 
motion. For the purposes of this paper, a system is considered stable if 
its motion is either damped or purely sinusoidal; only timewise divergent 
motion is considered to be unstable. In order to determine the stability 
boundary itself, attention can be restricted to simple harmonic motion. 
However, such simple -harmonic-motion analyses often yield a multiplicity 
of boundaries, and it is necessary to derive equilibrium conditions for 
^ divergent oscillatory motions in order to determine the degree of insta- 

bility in regions separated by the various boundaries and, thus, to 
identify the primary stability boundary. 




operator. 


V-4 inverse of 

Subscripts : 
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In line with the foregoing discussion, the method of approach is 
first to derive the equilibrium conditions for sinusoidal motion. These 
conditions are then extended to apply to divergent motion by means of 
analytic continuation. A Nyquist diagram technique is used to determine 
states of stability and stability boundaries . 


UNSTIFFENED CYLINDER 
Derivation of Equations 


It is assumed herein that the cylinder wall may deform into any 
number of sinusoidal waves around its circumference and into sinusoidal 
waves of any wave length and constant amplitude along Its length and, 
further, that the motion is simple harmonic in time; spacewise divergent 
motion (motion increasing in amplitude along the cylinder) is specifi- 
cally excluded. Thus, the lateral deflection of the cylinder wall may 
be written 


w(x,0,t) = Re^we“^ vx e^ a) ^cos nsj 

r -iv(x-! 

I T .TCI ' 


= Re 


we 

U 


cos nfi 


(n « 2, 3, 4 , . . .) 


CD 


where w" is the complex amplitude of the potion, v is the real wave num- 
ber of the longitudinal waves, n is the number of full waves around the 
circumference, and cs is the frequency of oscillation. This assumed 
deflection shape will be the basis for the determination of the air 
forces exerted on the cylinder, the equilibrium condition, and, through 
it, the criteria for flutter. Note that n = 0 and n = 1 have been 
specifically excluded from consideration in this panel-flutter analysis 
because neither of these two motions involves panel action. The first 
value n = 0 represents pure dilation or contraction of the cross sec- 
tion; the second value n = 1 merely represents a rigid-body translation 
of the cross section. 


Air forces .- The air forces exerted on an infinite cylinder vibrating 
harmonically in still air have been reported in the literature. (See, for 
example, ref. 13.) Although some of these .results can be applied to the 
present problem of determining the forces exerted on an oscillating cylin- 
der in moving air, it is convenient to derive this result directly. The 
unsteady-flow problem can be reduced to a steady-flow problem by means of 
a moving coordinate system in a manner similar to that of reference 7* 
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Accordingly, the perturbation pressures exerted by a steady external 
or internal flow of Mach number M]_ on a stationary cylinder deformed in 

the shape 

(_ -ivx-, \ 

w = Reive cos n9j (2) 


have been determined in appendix A. The assumed deformation (eq. (l) ) 
has the form of a wave traveling in the positive x-direction with veloc- 
ity cd/ v; hence, the flow of Mach number M outside the vibrating cylin- 
der is equivalent to a flow of Mach number M - au/vc outside the 
stationary cylinder plus a flow of Mach number -ca/ vCj_ inside the sta- 
tionary cylinder. Thus, by making the appropriate substitutions, the 
steady-flow results of appendix A are readily combined to give, for the 
outward air force exerted on a vibrating cylinder, 


z(x,0,t) = Rv^Re 


-iv{x- 22t\ 

Zwe v. v / cos n0 


( 3 ) 


« 


where 


l = 


-pc 2 F(M-k, v,n) 


+ 



(*) 


In equation (it-), 


F(|,b,n) 


I 2 

I 2 


KnCU 

£VU) 


for b| > 0 
for b£ < 0 


(Ul< 1)' 

> 

(III >i). 


( 5 ) 


r 




*i(|,b,n) 


2 

SVC5) 


1 2 ^n( z ) 

zJ n ’(z) 


(U I < 1) 


> 


(III > 1)J 


(6) 
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with. 


5 = M ~ a 2 

z = |h| J |2 - 1 


In equations (5) and (6), J n is the Bessel function of the first kind. 

Is the Bankel function of the first or second kind, and 3 ^ and 

are modified Bessel functions of the first and second kind, respectively. 
(See, for example, refs. 14 and 15.) The primed quantity in the denomina- 
tor of each function is the derivative with respect to the entire argument. 
Note that a dimensionless frequency k = cju/vc and a dimensionless wave 
number v = Rv have been introduced in equation ( 4 ). 

Equilibrium condition .- Donnell's equation for the equilibrium of 
thin, cylindrical shells (ref. 10} may be written in the modified form 
(see, for example, ref. 16) 


y^w + _4- J 1 2 1 vV 

R 


3Z ' " xxxx 


- h |^x v xx + ~ w ee) + Ps hw t.t - 1 (T) 


where the subscripts x, 8, and t on w indicate differentiation with 
respect to these variables. The operators and are defined by 


v» - (£- + i JLY 8 

\dx 2 R 2 8e 2 / 

and 

=: 1 


Substitution of w and l from equations (l) and (5) into equation (7) 
gives 


P + if + 3T 


12(1 - n 2 ) 

(I) 


v2 


+ h 


( f ° * + i ° 8 ) 


p s ha? - Rv 2 ! = 0 




y 
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as the condition for equilibrium of the motion of the cylinder vail in 
the presence of the moving air outside and the stationary fluid inside 
the cylinder. Written in terms of pertinent nondimens ional parameters, 
this expression becomes 


L = ?T |( k2 - A 2 ) + k ^ n ). - F (M-k,P,n) ‘ i |f|- €k * 0 


( 8 ) 


vhere 


c s h 

fi = — ■ ====== /v 2 (l + 

l/l2(l ~ 2) 



12(1 - u 2 ) 


1 . 

+ 5 + s£ 

(f 


-sr 

+ ° x + 3? ° 9 

a/ E, and 

Cs 

= /e/p s 

is the speed 


of sound in the cylinder material. Note that, in equation (8), an addi- 
tional term -i — ek has been added in order to include qualitatively 

Iv| 

the effects of structural damping of the Sezava (viscous) type. The 
damping coefficient e is actually related to the parameters of the 
system in a rather complicated fashion but is always positive. A pre- 
cise definition of e is not necessary because, in the following deriva- 
tion, it is considered to be a vanishingly small, positive quantity. 


Extension of Equilibrium Equation to Divergent Motion 

The flutter equation (eq. (8)) has been derived for sinusoidal motion. 
As has been pointed out, an extension to divergent motion (of the form of 
eq. (l) with o> a complex number having a negative imaginary part) is 
necessary. This requires the analytic continuation of equation (8) into 
the complex k-plane. The only terms in equation (8) for which this con- 
tinuation is not trivial are the air -force functions F and F-j_ . The 

analytic continuations of these functions are presented and discussed in 
appendix B. These functions are analytic throughout the half -plane corre- 
sponding to divergent motion and approach the values given by equations ( 5 ) 
and (6) as k becomes real. That the analytic extension of equation (8) 
into the lower half of the m-plane is indeed the equilibrium equation for 
divergent motion of the cylinder can be rigorously shown by the application 
of Fourier transform analysis. 
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Stability Determination by Use of Nyquist Diagrams 

In order to determine the stability states of the structure whose * 

behavior is governed by equation (8), it is expedient to examine the 
variation of the function L in the L-plahe as k traverses, in the 9 

clockwise direction, a curve enclosing that portion of the k-plane 
corresponding to divergent motion of the cylinder (the lower half -plane 
for positive values of v or the upper half -plane for negative 
values of v). This, of course, is the well-known Nyquist diagram. 

(See, for example, ref.' 7 .) The number of resulting clockwise encircle- 
ments of the origin in the L-plane is equal to the number of zeros of L 
minus the number of poles of L enclosed by the curve traversed in the 
k-plane. Poles of L do not occur in the unstable half -plane since the 
functions comprising L are analytic everywhere in that region, (it 
should be noted -that, in all instances, the real k-axis is excluded from 
the unstable half -plane.) Therefore, the Nyquist diagram gives directly 
the number of roots of equation (8) for which the frequency of oscillation 
has a negative imaginary part; in other words, it gives directly the num- 
ber of modes of unstable motion. 

With regard to the practical application of the Nyquist diagrams in A 

the present case, it should be pointed out that actual computation of 
successive values of L is ■unnecessary; that is, through careful examina- 
tion of the nature of the functions comprising L, 'it is possible to con- 9 

struct approximate diagrams correct in all essential features and to see 
clearly the conditions -under which critical changes (corresponding to 
flutter boundaries ) occur. Further, it is only necessary to consider 
positive values of v because, if v is replaced by its negative and 
k by its conjugate, then L is replaced by its conjugate; therefore, no 
new flutter boundaries would result from the negative values of v. 


Analysis of Empty Cylinder 

Flutter criteria .- If the cylinder is assumed to be empty (p^ = 0), 
equation (8) reduces to 

L s £S- |(k 2 - fi 2 ) - F(M-k,v,n) - i y^y ek = 0 (9) 


Let the damping coefficient e be very small but positive. With this 
restriction, the functions comprising L T (for v > 0) vary with k in 
the manner shown in figure 2. Note that the inclusion of the infinitesi- 
mal damping is influential only in the range M-l<k<M+l where 
Im(-F) = 0; elsewhere, Im(-F) is large in comparison with the damping 
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term. In the range M-l<k<M+l, the damping makes Im(L) nega- 
tive for positive k and positive for negative_ k. Hence, for subsonic 
flow (fig. 2(a)), the damping insures that Im(L)_ = 0 only at k = 0; 
for supersonic flow (fig. 2(b)), it permits Im(L) = 0 only at k = M - 1 
(in the limit as e — *0) . The significance of this will be made apparent 
by considering the variation of L as k traverses the path shown in 
figure 3 (with the other parameters held constant). 

First, consider subsonic flow. (See fig. 2(a).) Two possible 
resulting variations of L are shown in the following Nyquist diagram 
(where the full circle at infinity corresponds to the infini te semicircle 
in the k-plane ) : 



For different values of the parameters, paths similar to I or II may 
be traced. For path I, there is no encirclement of the origin and the 
cylinder is stable; for path II, oh the other hand, one encirclement 
occurs and the cylinder is unstable. Since, by virtue of the damping 
term, 3m(L) must pass through zero at k = 0, it is apparent that the 
boundary between these two conditions is a static (divergence) bo undar y 
defined throughout the subsonic range by 


Re (L) k=0 = 0 = -^- | fi 2 - F(M,v,n) 


W 


( 10 ) 
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Equation. (10) may be put into the form 



( 11 ) 


where 


12(1 - y.2) 


1 

n 2 

v 2 ! 

'1 + 

\ V 2 / 

-V 2 

i + 

V 

2 + °x + zs °e 


and. 


B = -Lgil ~ H 2 ) g(M,v,n) 

wf 52 


eire positive numbers. Equation (ll) has the one real root 



Thus, for selected values of the other parameters, the thickness ratio 
corresponding to the stability boundary in the subsonic range may be com- 
puted directly. This, is the bounding value of h/R above which the 
cylinder is stable. For certain ranges of the parameters, namely when 
the Mach number is only slightly less than 1 or when the wave length of 

flutter is very large, the approximation F(M,v,n) is permissible. 

When this approximation .is valid, equation (10) may be solved directly 
for M to yield 




v 2 ( 1 - M 2 ) 
2n(n - 1) 


« 1 


( 13 ) 


Equation ( 13 ) gives the approximate bounding value of M below which 
the cylinder is stable. 
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For supersonic flow (fig. 2(b)), tbe possible Nyquist diagrams 
resulting when k traverses the path of figure 3 appear as shown in the 
following sketch: 



Since, now, the sign of the imaginary part of L always changes at 
k = M - 1, the criterion for a boundary is seen to be 


Ee ( L )k=*l-l = 0 = y- |[(M - l) 2 - ft 2 ] - F(l,v,n) (14) 

Thus, the instability is dynamic, and, since k = the flutter mode 

is a traveling wave whose propagation velocity ^ is the velocity of 
external flow minus the velocity of sound. 

The solution of equation (14) for h/n can be carried out to yield 
a result in the same form as equation (12) with appropriate redefinition 
of A and B. In the supersonic case, however, the F-f unction is inde- 
pendent of M and is actually equal to -l/n; hence, equation (l4) may 
be solved directly for the Mach number of flutter to give the convenient 
form . 
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Determination of critical stability boundary .- In computing sta- 
bility boundaries, a factor which must be kept in mind is the admissi- 
bility of all longitudinal wave lengths and of all integral values of n 
greater than or equal to 2. The critical, boundary must necessarily 
correspond to the values of v and n for which the cylinder is most 
prone to flutter. This critical boundary can be obtained either by 
expressing the bounding value of h/R as a function of the rest of the 
parameters (as in eq. (12)) and maximizing this value of h/R with 
respect to both v and n or by finding the bounding value of M (as 
in eqs. (13) and (15)) and minimizing with respect to v and n. Of 
the two alternatives, the latter is more easily accomplished, analytical 
minimi zation of M being possible. For subsonic speeds, however, the 
expression for M (eq. (13)) is only approximate and, in some instances, 
it becomes necessary to maximize h/R as given by equation (12). This 
is most readily achieved by graphical methods . 


Minimization of M can be performed by first minimizing fi 2 with 
respect to v and then minimizing the resulting M with respect to n. 
The minimization with respect to v is particularly simple when there 
are no imposed midplane stresses = Cq = 0); it can be performed with 

respect to the quantity . The result is 


at 


^min 2 ~ 



/l2(l - p) 2 


(16) 


( v + ff . & . (17) 

R 

However, the quantity I s Itself a minimum at v = n; hence, 

equation (17) can be satisfied only when 


^- 2 (.j- ~ > ixn 2 (18) 

R 


and, if this condition is violated, a 2 can never achieve the value 
given by equation (l6). The minimum value it can achieve is 




12(1 - p 2 ) 


+ 



(19) 
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If the result given by either equation (l 6 ) or equation ( 19 ) is sub- 
stituted into equation (15), the minimum value of M is seen to corre- 
spond to the smallest admissible value of n (that is* n = 2 ); hence, 
for supersonic flow, the critical boundary is given by 



Similarly, for subsonic flow, if the approximate equation, equation ( 13 ), 
is valid, the substitution of equation (16) into equation (13) yields the 
critical boundary 



Numerical example .- For Illustrative purposes, the critical stability 
boundary (n = 2) has been computed for an empty, unstressed aluminum cyl- 
inder at sea level ; additional curves corresponding to n = 3> and 

10 have been obtained for comparison. These results are shown by the 
solid curves in figure 4 . Portions of the curves in the subsonic range 
were obtained by graphically maximizing h/R with respect to v as pre- 
viously mentioned. It is interesting to examine the wave lengths of the 
flutter modes corresponding to these boundaries. In the range where equa- 
tion (l8) is satisfied, solution of equation (17) yields two different 
wave lengths for the same critical thickness ratio. The larger one of 
these two wave lengths is shown in figure 5 for supersonic Mach numbers 
in the form of a plot of the aspect ratio of the flutter mode n/v (the 
ratio of longitudinal to circumferential wave length) against Mach number 
for n = 2 , 3 , 4 , and 5. The aspect ratios associated with the smaller 
wave lengths are not plotted; they are merely the reciprocals of the ones 
shown. At the higher Mach numbers, the critical value of h/R becomes 
large enough to cause equation (l8) to be violated; when this happens, 
the two flutter modes coalesce to give a single flutter aspect ratio of 
unity as shown by the horizontal cutoff line in figure 5* 
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Ab has been pointed out, Donnell's theory is somewhat inaccurate for 
small values of n. (See, for example, ref. 12.) In order to obtain some 
idea of the magnitude of the resulting error j a new ^-function based on a 
simplified version of Flugge's cylindrical-shell theory (see, for example, 
refs. 11 and 12) has been derived and minimized with respect to v; the 
details of this derivation are presented in appendix C. ■ Computations in 
the supersonic range with the new (Fliigge) fi-f unction result in the dashed 
curves shown in figures k and 5. As in Donnell's theory, the Fliigge 
theory yields, for the lower thickness ratios, two values of wave length 
for which Q. (and, hence, the Mach number) is a minimum. In contrast 
with Donnell's theory, however, the two wave lengths are not equally crit- 
ical ; the larger one always yields a lower Mach number. For this reason, 
only the higher aspect ratio and its associated stability boundary are 
plotted in the figures. The critical boundary (n = 2) is found to require 
thickness ratios approximately JO percent higher than those predicted by 
Donnell's theory. The two curves for n = 3 still differ by 10 to 15 per- 
cent; but, for n = k, the two theories agree very well. For the sake of 
clarity, the stability boundary associated with the lower aspect ratio has 
not been shown in figure Ij-. It should be remarked, however, that this 
boundary agrees very well with the boundary given by Donnell's theory, 
even for n = 2. The practical implications of the results shown in fig- 
ures 4 and 5 are discussed subsequently. 


Analysis of Fluid-FillecL Cylinder 


When the cylinder is assumed to contain a fluid, the additional term 
p i ('cf’) *i^r" in equation (8) must be included. A typical plot of 


this function for real values of k is shown in figure 6. Note that the 
force is always real and becomes infinite at the resonant frequencies 


kj = irJ 1 + (j - 1, 2, 3, . . 0 (22) 


where the Zj's are the zeros of J n (z). 

Since Fj_ is always real for real values of k, the imaginary part 

of L can again be zero only at k = 0 In the subsonic range and at 
k = M - 1 in the supersonic range. Stability boundaries can consequently 
occur only at these frequencies. In the subsonic range, the inclusion of 
Fj_ has no effect upon the stability boundaries since F.^ = 0 at k = 0. 
In the supersonic range, however, the inclusion of Fj_ has a considerable 
effect on the stability boundaries, not only because of the additional 
force at k = M - 1 but also because of the resonances of F-j_. 
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The condition for the stability boundary at supersonic speeds is 

Ke ( L )k=M-i - 0 - ~ |[( M - i) 2 - n 2 ] + i + (23) 

This equation can be solved for h/R in the same form as equation (l2)j 
that is , 


1 + 


HE aA 1/5 j .(b_ IW~lE \ 1/3 


R 2 P 27 . 


+ — - *7- + — - 

2 J If- 27 


(24) 


where now 


and 


B = 


. 12(1 - n2) 

1 

v 2 |l + 

Wl + 

A 2 

[ V 2 / 

L v 

v2/ 

12(1 - Ji2) 

1 

1+ f 


n 2 = 


+ c x + - 


(M - l) 2 


m- 


t 1 + ST 


V)*" 




The variation of h/R with Mach number for particular values of 
the other parameters is shown in the following sketch: 
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Note that h/R hecomes infinite when the Mach number is equal to 1 + kj 
(j = 1, 2, 3> . . .)• (See eq. (22).} This arises from the infinite- 
ness of Fj_ and, consequently, of B at these Mach numbers. Note also 
that the footpoints of the secondary stability boundaries occur when 
B = 0. The numbers within the regions separated by the solid lines indi- 
cate the degree of instability as determined by use of Nyquist diagrams. 

Since the value of h/R for the empty cylinder is given by equa- 
tion (24) with F-j_ equal to zero, and since Fj_ is positive for values 

of M less than 1 + k]_, it can be seen that the fluid inside the cyl- 
inder has a destabilizing effect. 

In order to find the critical stability boundaries, all values of 
v and n must be considered. Of particular interest in this connec- 
tion is the fact that, for very large values of v, the infinities shown 

/ 

in the sketch all approach a value of M equal to 1 + (See 

eq. (22).) Thus, the critical stability boundaries would appear as 
shown in the following sketch: 


h 

R 



C4 

For Mach numbers greater than 1 + the cylinder is unstable to an 

infinite degree. If the fluid inside the cylinder is air at the same 
temperature as the surrounding air, this limiting Mach number would be 
equal to 2. If the cylinder contains a relatively incompressible fluid, 
however, this Mach number could be very high. In any event, the result 
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is somewhat anomalous j it is probably caused by the use of linearized 
potential -flow theory and, undoubtedly, would not occur for real fluids. 


Seme Remarks About the Solution for 
an Unstiffened Cylinder 

The stability criteria for the infinitely long, unstiffened cylinder 
which were derived in the preceding sections were obtained by including 
the effects of structural damping and then taking the limit as the damping 
approached zero (e — ^0) . This procedure was followed because different 
criteria are obtained when damping is not considered (e = 0). This impor- 
tant fact is illustrated in appendix D where the stability criteria for 
the empty cylinder with zero damping are discussed. Since structures 
always exhibit some damping, it is apparent that stability criteria 
obtained by taking the l imi t as e approaches zero are more realistic 
than those obtained by setting e identically equal to zero. It is 
interesting to note that, in this case, the addition of damping makes 
the structure more prone to flutter. This somewhat surprising result 
may be explained by the fact that a damping force, even though in itself 
dissipative, can cause phase changes in such a manner as to allow the 
moving outside air to feed more energy into the structure; the result is 
a net energy gain. 

With regard to the applicability of the results for the infinitely 
long, unstiffened cylinder to a flutter analysis of cylinders of finite 
length, the following remarks are in order. 

It Is clear that the results for the infinitely long cylinder would 
be applicable to a finite cylinder only if the wave lengths of the flutter 
modes were small in comparison with the length of the finite structure. 
But, the most critical wave lengths for the infinitely long, unstlffened 
cylinder are very large. (See fig. 5-) It is conceivable that, for a 
finite structure, the flutter mode would tend to settle on the smaller of 
the two possible longitudinal wave lengths discussed previously. However, 
it should be mentioned that, for higher Mach numbers (above M = 5)> even 
these smaller wave lengths are fairly high; for n = 2, for example, wave 
lengths from one to three times the radius would be experienced. 

Although the analysis has been carried out for the case of an 
unstiffened cylinder, the flutter criteria may also be applied to the 
case of a cylinder with essentially rigid, longitudinal stiffeners. 

These longitudinal stiffeners would have the effect of raising the mini- 
mum value of n and, hence, of decreasing the critical thickness ratio 
of flutter. For example, the curve for 20 stringers (n = 10) has been 
plotted in figure b. Also, the decrease in circumferential wave length 



20 


NACA TO 3658 


would cause an attendant decrease in the longitudinal, wave length at 
flutter and so might possibly make the analysis applicable to practical, 
finite cylinders. 

The effects of midplane tension stresses (caused, say, by a static- 
pressure differential across the cylinder wall) have not been investi- 
gated quantitatively; it is evident, however, that these stresses would 
increase the critical Mach number since they always increase ft (defined 
immediately after eq. (8)). 

In addition, a qualitative investigation of the behavior of the 
ft-function for Donnell's theory shows that the longitudinal stress a x 

has no effect on the critical wave length but that the circumferential 
tension Oq tends to decrease both the wave lengths corresponding to 

minimums of ft. More importantly, the hoop tension makes the lower of 
these two wave lengths critical because it is less influential for the 
lower wave lengths. Thus, for large circumferential stresses, in view 
of the aforementioned agreement between Flugge's and Donnell's theories 
for the lower wave lengths, the preceding analysis, based on Donnell's 
theory, might indeed give useful results for finite cylinders* 

This completes the discussion of unstiffened cylinders; an analysis 
of ring-stiffened cylinders is presented in the next section. 


RING-STIFFENED CYLINDER 


The ring-stiffened cylinder consists of the unstiffened cylinder 
with added, rigid ring stiffeners which prevent radial deflection at the 
locations x = ±ja (j = 0, 1, 2, . . .). (See fig. 1(b).) The stiff- 
eners are assumed not to interfere with the flow of air outside or of 
fluid inside the cylinder. The analysis proceeds along the lines of that 
in reference 7* 


Derivation of Equations 

If the assumptions are made that the cylinder wall may deform into 
any number of sinusoidal waves around its circumference and Into any shape 
periodic over 2 bay lengths in the longitudinal direction, the deformation 
may be written, with complete generality, as 


w(x,9,t) 



— l- 




mjtx 
a Icct 

e cos n0 


( 25 ) 
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provided that the coefficients 


aja satisfy the constraining relations 



09 

H 

m=-co 


► 

(-1)% - 0 


(26) 


These conditions correspond to zero deflection at the ring locations . 
They may also he written as 



m=~oo 


a a “ 


(m odd) 

CO 


0 


\ 





(m even) 


(27) 


Only one circumferential term is included in equation (25) because there 
is no structural or aerodynamic coupling between the various cosine terms. 
Also, in this case, n = 1 is admitted. The assumption of periodicity 
over 2 bay lengths is made because it permits a considerable degree of 
generality without overcomplicating the analysis. It is believed that 
the critical flutter mode would be of this type. 

Air forces . - In view of the linearity of the aerodynamic theory, 
the air forces exerted on a cylinder executing the motion given by equa- 
tion (25) may be determined by separately considering each term of the 
summation over m and superposing the results. Hence, the aerodynamic 
loading on a cylinder deforming in the shape of the traveling wave 


w m.( x , 0 >t) = Re 



cos nS 


(28) 


is sought. 
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If equation (28) is compared with equation (l), it becomes apparent 
that the air force is given by equations (3) to (6) when the substitu- 
tions w = a m and v = mr/a are made. Thus, the total air force may 

be written immediately as 


z(x,e,t) 



■i max 

, 1 eT" itDt 

Zjjja^e e cos n0 


where, for m ^ 0, 



PiCf^Fi^ |,ma,nj 


( 29 ) 


( 50 ) 


The functions F and F-^ are given by equations (5) and (6). They are 
repeated here for convenience. 


F(|,b,n) = £ 2 


^'(C) 


(Ul < !)1 


H^J)(z) 

( , _ fl for bS > 0] 

zn^y(z) 

\ ]2 for b| < 0) 


> ( 51 ) 


(lei >DJ 


F±(S >b,n) 


2 

£V(£> 


I 2 


J n (z) 

ZJ n ’(z) 


where 


(III < 1)1 


> ( 32 ) 


(ill >Dj 


i = M - 1 2 

z = |b J sji 2 - 1 
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In equation (30), k = <m/itc and a = Rit/a-. The quantity a is the 
ratio of one-half the circumference to the distance between the rings. 
For the special case of m = 0, a limiting process gives 


l 0 = -pc 2 G(k,a,n) + k,a,n) 


(33) 


In equation (33), 


G(£,a,n) 


2 H^ J) (a|tl) 
aWH^'CaUI) 


11 for £ < 0 

12 for | > 0, 


(3*0 


Q^(S n ) 


= g2 J n( a U I ) 
a U! J n , ( a UI ) 


(35) 


Equilibrium conditions . - By virtue of the assumed longitudinal 
spatial periodicity of the deformation, satisfaction of equilibrium 
over any interval of length 2a assures satisfaction of equilibrium 
over the whole length of the cylinder. Thus, it is sufficient to write 
Donnell's equation for the segment 0 ^ x< 2a as 


v\r + -A- -*-2(3- - H 2 ) V7-4, 

R^ (h\ 2 

VrJ 


w. 


xxxx 


- 11 Px w xx + p- w eeJ + 


Ps^tt = J (x,e,t) + 2 q( 9 , t)8(x) + P 1 (0,t)5(x-a) 


(36) 


where Pq and P-]_ are the reaction forces exerted on the cylinder wall 

by the two ring stiffeners included in the interval and where 5(x) is 
the Dirac delta function. 
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P Q (e,t) = Re ^PQe^cos n 9^ 
P 1 ( 6 ,t) = Re^e^coe n0^ 


If P 0 and from equations (37) , v from equation ( 25 ), and l from 

equation ( 29 ) are substituted into the equilibrium equation, equation ( 3 6 ), 


and if this equation is then multiplied by e 
interval, there results 


and integrated over the 


2 a s D 


i 2 * 2 , n 2 \ 2 . 12(1 - H 2 ) 


/m.2^2 , n2 \ 2 


+ h^£ |_<j x + 2 §<r 9 


p s ha^ - B *=. jL % = P 0 + (-l) m P L 


Hence, for 0 m £ 0, 


P Q + (-l) m P X 

2n 2 c 2 pR 0 

£ 


where, for m ■£ 0 , 


(m = 0, ±1, ±2, ±3 ; • • ■) (58) 


P = ZB.UI 

p m p Hi 


(% 2 - k 2 ) - m 2 Fi^ |,ma,n) + m 2 F^-|,ma,n) + i jfj- e m k 
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with 


— 


. g - B . 


^ 12(1 - p. 2 ) 


• + -±±» 
\n maj 


12(1 - u 2 ) 

1 

OJ 

n 2fact + J2J1 2 

\ n mo,/ 


+ a x + 


\ma ) 


2 _ 

^0 


Note that, as in the case of the unstiffened cylinder, a term has been 
added to each of the J^'s to approximate the effect of damping (e m > 
In the limit as m — * 0 , 


g 

0 o = !( n 0 2 - k2 ) - ^(if) °i^ k,a,n^ + G(k,a,n) + ie^ 


where 


Oo - 


c s h 
£-fi 


V12C1 - v- 2 ) 


a u 2 ^ 12(1 - n 2 ) - 


a,. 


[n. + 


. . r 

\ R / 


Finally, substituting from equation ( 39 ) into the constraining relations 
equations (27), gives the following conditions for the existence of the 
motion defined by equation ( 25 ): 


( F o - ?x) 


m=-co 

(m odd) 


~ = 0 
0 m 


(to) 


( p ° + p i) £ r-° 


OH) 


111= -CO HI 

(m even) 


Inspection of these relations indicates that simultaneous satisfac- 
tion of the equilibrium equation and the restraint conditions can be 
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achieved in a nontrivial fashion in several independent ways. One possi- 
bility is that 



(which implies that = 0 when m is odd) and 



(m even) 


and another is that 


(te) 



(which implies that a^j = 0 when m is even) and 


00 



m=-oo 



(m odd) 


0 


(^ 3 ) 


As has been noted, the derivation herein and the results achieved 
exactly parallel those obtained in reference 7 for the infinite, flat 
plate. As pointed out therein, still another nontrivial solution is 
found when the restriction ^ 0 (which characterized the foregoing 
results) is removed. Specifically, if two or more j^'s of the same 
type (m odd or even) vanish simultaneously, a flutter mode may exist. 

Equations (42) and (43) are conditions for the existence of motion 
of constant amplitude; equation (42) corresponds to motion which is 
identical in each bay, whereas equation (43) corresponds to motion having 
the same amplitude from bay to bay but with alternating direction. 


Stability Boundaries 

As in the case of the unstiffened cylinder, the Nyquist criterion 
can be used to investigate the states of stability of the stiffened 
cylinder. In this section are given the results of a limited investi- 
gation which was chiefly concerned with the examination of the 
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following three-term approximation to equation (42): 


N = — + — + — = 0 

Pq ^2 P-2 


(44) 


Subsonic flew .- In applying the Nyquist diagram technique to the 
function N, account must be taken of the poles of N at the zeros of 
p 2 and p_ 2 - (There are no zeros of (3q.) The results of such an 

application with e — >0 show" clearly that, for subsonic Mach numbers, 
the only possible instability is divergence; further, the three -term 
approximation leads to as many as three static stability boundaries as 
shown in the following sketch: 



M 


In this sketch, a typical variation of (h/R) max with M (with other 

pertinent parameters fixed) is shown. The subscript max indicates that 
the thickness ratio is maximized with respect to n. The numbers indi- 
cate the n umb er of unstable roots of equation (44) corresponding to each 
region. Note that the upper curve corresponds to the condition that two 
or more f^’s of the same type vanish simultaneously, whereas the others 

curves result from the condition (nJ^-q = 0* 

Additional curves would result from the addition of more terms to 
N and still more curves would result from the consideration of the odd 
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solution, equation ( 43 ); further, examination of a four-term approxima- 
tion (m = 1, -1, 3, and -3) to equation ( 4 - 3 ) indicates that, for reason- 
able values of a, the most critical boundary for either the even or the 
odd solution must stem from the condition = l)k=0 = (^ or 


very s mall values of cx, the distance between rings becomes large and 
values of m other than -unity may be critical.) The position of this 
critical boundary has been computed for an unstressed aluminum cylinder 
at sea level for various representative values of the aspect-ratio 
parameter a. The results are shown in figure 7 - These curves corre- 
spond to various values of n ranging from 5 upward as illustrated on 
the plot; since n is large, Donnell’s theory is sufficiently accurate. 


Note that the thickness ratios required for static stability at 
subsonic Mach numbers are extremely small; therefore, divergence at 
subsonic speeds is probably not a critical design factor. 


Supersonic flow . - Application of the Nyquist criteria for supersonic 
flow is not so readily accomplished without the performance of further 
computations. It is clear, however, that divergence boundaries extend 
into the supersonic range and that flutter boundaries arise which probably 
become more critical than the divergence boundaries. 


The def ini tion of these flutter boundaries requires that solutions 
be obtained to equations ( 4 - 2 ) and ( 4 - 3 ) or to suitable approximations, 
such as equation ( 44 -). This, in turn, requires extensive computations, 
especially since the resulting thickness ratios must be maximized with 
respect to n. This maximization can probably only be achieved laboriously 
by making each computation for several values of n. 


CONCLUDING REMARKS 


A preliminary theoretical investigation of the aeroelastic stability 
of infinitely long, thin-walled unstiffened and ring-stiffened circular 
cylinders has been conducted by using Donnell's cylinder theory and 
linearized unsteady potential -flow theory. A limited study of the 
resulting stability criteria has yielded the following information. 

For unstiffened cylinders with vanishingly small structural damping, 
the only possible instability at subsonic Mach numbers is static diver- 
gence; in supersonic flow, however, flutter is found to occur for suffi- 
ciently thin cylinders in the form of a traveling wave whose propagation 
velocity is the velocity of the external flow minus the speed of sound. 
For an empty, unstressed aluminum cylinder at sea level, the critical 
boundary is found to correspond to a mode of deformation having only two 
waves around the circumference. For this case, the use of the more 
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complicated cylinder theory of Fliigge may be necessary. Relatively 
large ratios of cylinder -wall thickness to radius are found necessary 
for stability, and the wave lengths corresponding to flutter are found 
to be very large. With the addition of large, circumferential mid- 
plane tension stresses, however, the thickness ratios and wave lengths 
would be reduced. 

The addition of an internal fluid has a destabilizing effect on the 
unstiffened cylinder. In fact, the anomalous result is found that for 
Mach numbers greater than 1 plus the ratio of the speeds of sound in the 
fluids inside and outside the cylinder, no adjustment of the physical 
properties of the cylinder will render it stable. 

The presence of even the smallest amount of structural damping is 
found to be an important factor in analyses of infinitely long, unstiff- 
ened cylinders . 

For the ring-stiffened cylinder it is found that flutter is not possi- 
ble at subsonic Mach numbers and that only very small thickness ratios are 
required to prevent divergence. Although both panel flutter and divergence 
are possible at supersonic Mach numbers, no numerical results have been 
obtained; extensive computations would be required for a complete determi- 
nation of the stability boundaries in this range. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 11 , 19 56 . 
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APPENDIX A 

STEADY LINEARIZED PLOW PAST A STATIONARY DEFORMED CYLINDER 


In this appendix axe derived the perturbation pressures exerted by 
a steady external or internal flow of Mach number on an infinitely 
long thin-walled cylinder with the deformation 


w 


Re 




(Al) 


where n is a positive integer and the wave number v may be either 
positive or negative. The Mach number Mq is positive for flow in the 

positive xq-direction and negative for flow in the negative xq-direction. 
(See fig. 8.) 

In the cylindrical coordinate system (xq,r, 9 ) , the linearized 
potential-flow equation is 


(l - Mq^)$xqXq + 0 rr 


+ ? + “ 0ee - 


0 


(A2) 


where <f> is the velocity potential and the subscripts xq, r, and 8 

indicate differentiation with respect to these parameters. With the 
assumption that 


0(xi,r,8) 



(A3) 


equation (A2) becomes 


f rr + ~ f r + 



Uk) 
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which, for supersonic flow (|Mij >1), has the solution. 

f(r) = mi 1] (|v| ^M X 2 - 1 rj + mi 2) ^|v| ^Hq 2 - 1 
and, for subsonic flow (|Mi| < l), has the solution 

f (r) = AI n ^|v| \jl - rj + BK n (jv ( \j± - M X 2 rj 


Air Plow Outside Cylinder 


(A5) 


(A6) 


For air flowing over the outside of the cylinder, the formulation 
of the problem is completed by the specification of the boundary condition 


$ r Ui,R,e) = M;]CW Xi 


(AT) 


at the cylinder wall and the proper conditions at infinity. 

By use of equations (Al) and (A3), equation (A7) becomes 


f r (R) = -lM lC vw 


(A8) 


The resulting pressure perturbation Zp can be calculated from the 
linearized Bernoulli's equation 


Zp(x L ,r,0) = -pM 1 c^ Xi 


(A9) 


and is 


^4 ) ( x 1 ,r,9) 


pMqcRe^ifve ^ VX -'-cos nsj 


(A10) 
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Supersonic flow .- For supersonic flow, the Sonmerfeld condition 
(that is, that there be no incoming disturbances from infinity) requires 
that, for large values of r, the velocity potential $ be essentially 
a function of x^_ - pr for > 0 and a function of x-j_ + pr for 

Mp < 0 ^ where p = \/mi 2 - l). Then, from the relation for (eq. (A3)), 
it is apparent that, for large values of r, f(r) must behave like 
Mi 

i 'Rf vr 

e . Substituting the asymptotic approximations (ref. 14) 




~ ji 

y rtz 


Hn '(*> ~ 


■-Hk 


i -i[z-^n+i)|J 


(All) 


into equation (A5) gives, for large values of r. 


OZ e - 1 ( n+ i)l A e i l v|pr + ^ 

it Ivipr [_ 


(A12) 


Hence, the Sonmerfeld condition requires that B = 0 for Mqv > 0 
and A = 0 for Mqv < 0. Therefore, equation (A5) becomes 


Ivh/Mi^ - 1 r 


1 for Mqv > 0 

2 for Mpv < 0 


If, now, A is evaluated by use of the boundary condition, equa- 
tion (A8) , and the resulting expression for f is- substituted into 
equation (AlO) , the perturbation pressure on the cylinder wall is 
obtained as 
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Ap(xq,R, 8) = pc?Re 


r(j). 


Rv% x 2 


— ^ lcos EL0 


zHq (z) 


we 



1 for M]_v > 0 

2 for M^_v < 0 


(A14) 


where 


z 



Subsonic flow. - For subsonic flow over the outside of the cylinder, 
the potential must remain finite at infinity; hence, equation (a 6) reduces 
to 


f(r) = BK n 




(A15) 


If B is evaluated by using the boundary condition, equation (A8) , 
and the resulting expression for f is then substituted into equa- 
tion (A10), the resulting perturbation pressure on the cylinder wall is 


£p(x-j_,R,0) = pc%le 


Rv%i 2 


KnU) 

5V(5) 


— -iVXn 

we J -cos n0 


(Al6) 


where 


i = 



Fluid Flow Inside Cylinder 

It is assumed that the fluid inside the cylinder may be other than 
air and may exist under different conditions from the air outside the 
cylinder. Thus, the boundary condition at the cylinder wall is 
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f r (R) = -iM^c^vw 


(AIT) 


and Bernoulli's equation may be written as 


^3 i (x 1 ,r,0) 


Pi M l c i Re ( if ve” ivx lcos 



(A18) 


where the subscript i refers to the properties of the fluid inside the 
cylinder . 


For either supersonic or subsonic flow, the velocity potential must 
remain finite throughout the cylinder. If this requirement is to be met 
at r = 0, the solutions given by equations (A5) and (a 6) must reduce, 
respectively, to 


and 



- (A19) 


f (r) = Aln^i v | V 1 " M l 2 


(A20) 


With the constant A determined through the use of the boundary condi- 
tion, equation (A17), and the resulting expressions for f substituted 
into equation (Al8) , the perturbation pressure for |M-]_| > 1 becomes 






1?— L ve- lvx W n0 
zJ n '(z) 


and, for 


|“l| < 


f5> 1 (x 1 ,R,0) - p 1 c i d Re 


Rv%i 2 


In(0 - we- ivx lcos n9 

5V (5) 


(A21) 


(A22) 


These results hold, for flow in either direction. Note that, for these 
air forces, the pressure is always in phase with the motion. 
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APPENDIX B 

AIR FORCES ON A CYLINDER IN DIVERGENT OSCILLATION 


The crutvard lift force on a cylinder wall executing the simple 
harmonic motion (o> real) 


w(x,9,t) = Re^we -ivx e ldjt cos n oj 


(Bl) 


in the presence of moving air outside and a stationary fluid inside the 
cylinder has been determined. (See eqs. (3) to (6).) The amplitude of 
this force is expressed by 


l = -pc%^M-iL,Rv,n^ + Pi^i^Rv,*) 


(B2) 


In equation (B2) 


F(|,Rv,n) = | 2 


= I 2 


%(!) 

ivTET 

Hn J) (z) 

zEn^ (z) 


J = 


[1 for Rv£ > 0 
2 for Rv| < 0 


(Ul < 1) 


(III > 1) 


> (B3) 


and 


, . 2 ^(S) 

Fj (|,Rv,n) = 5 — — 

1 ’ lln’(l) 


(III < 1) 


r W 


= 1 2 


Jn(z) 

zJn* (z) 


(III > D 
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where 


6 = R|v|^ 1 - i 2 

z = R|v|^£ 2 - 1 

It is desired to extend these air-force functions to apply to 
divergent motion of the cylinder (a> complex with a negative imaginary 
part) . This extension can be made by using analytic continuation; in 
fact, it can be shown by Fourier transform analysis that the necessary 
and sufficient conditions for the extensibility of the air-force func- 
tions to divergent motion are that the air-force functions possess analytic 
continuations which have no singularities anywhere in the unstable half- 
plane (lm(u>) < 0) and reduce smoothly to equations (B3) and (b 4 ) on the 
real axis. Singularities may possibly exist along the real axis, but 
the real axis is specifically excluded from the unstable half-plane. 

Consider the function F(|,Rv,n). Examination of the manner in 
which appears in F (see eq.. (B 2 )) shows that the unstable half of 
the £ -plane is the upper half -plane if v is positive and the lower 
half -plane If v is negative. It can be verified that the desired 
analytic continuation of F into the proper h a l f -plane 1 b given as 
follows : 

I 2 k i1 (r|v|\/i - I 2 ) 

F(|,Rv,n) = ^ 7 f - — y (® 5 ) 

r|v| \Ji - i 2 k^vrm y 1 - 1 2 J 

where the desired branch of the multiple -valued function F is the one 
for which the cuts are as shown in the following sketch: 




v>0 


v<0 
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and 


F( fz, Rv,n) 


2H^ 1} (R|v|) 
RIvlH^ 15 *(R|v |) 


24 2) (R|v|) 

(p ) « 

R|v|H^ ; (R|v|) 


(v > 0) 


> (B6) 


(v < 0) 


As Im(|) approaches zero from, the proper direction, the definitions 
given by equations (B3) result. Furthermore, no branch points or cuts 
of the function F occur in the unstable half of the l-plane, and, in 
addition, Nyquist diagram techniques can be employed to show that there 
are no poles of F in this region. 


In a similar manner, the analytic continuation of F^ can be 
expressed in the form 


F ± (S,Hv,ii) 



R|vl\/l - |2 l n ^R| v l\Jl - i 2 j 


(B?) 


No difficulties with branch points occur for F-j.; the functions I n 
and In' are entire functions and the combination I n /z 1^ is an even 
function. Therefore, F^ is inherently a single-valued function of | . 

The air forces used in the analysis of the ring- stiffened cylinder 
can be similarly written for divergent motion. For m / 0, it is only- 
necessary to replace v by nur/a in equations (B5) to (B 7 ) • For m = 0, 
it is necessary to specify the air forces separately. Thus, for divergent 
motion (see eqs. (33) to (35 ) )> 


Iq = -pc 2 G(k,a,n) + PiCi^iJ— k,a,n] 

,c ± 


(B8) 
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where 


0(6, a, n) 


(at) 

azi 2) '(a|) 


(B 9 ) 


%(l,a,n) 


£J n (a6) 

l 

aJ n (a|) 


(BIO) 


The desired branch of the multiple-valued function G is the one for 
which the cut extends from the origin to infinity along the positive 
imaginary axis and 


G(l,n,a) 


4 2 ) U) 

aE^ ' (a) 


(Bll) 
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APPENDIX C 


EQUILIBRIUM OF UNSTIFFENED CYLINDER - FLUGGE* S EQUATION 


The simplifying assumptions used in deriving Donnell's differential 
equation (eq. ( 7 )) result in a loss of accuracy when the condition 
n 2 » 1 is violated; further, terms which become important when the 
long it udinal wave lengths become large (that is, when the behavior of 
a cylinder approaches that of a ring) have been omitted. 

The more complicated theory proposed by Flugge (see, for example, 
ref. 11) is not characterized by either of these limitations. Flugge* s 
equations may be written in the form of a single equation in w, which, 
for (h/R ) 2 « 1, reduces (see refs. 11 and 12) to 


D1 


vV + — 


12(1 - P 2 ) 
(h/R) 2 


w- 


xxxx 


2(2 - p) 

R 2 


1 

w xx 9 e + -r w eeee 
R^ 


+ 


^^xxxxxx + 6w xxxx00 + ■ " ---• W XX0000 + ~t w 000000 

R 2 R^ 


J 


J> + P B hw tt a 1 
(Cl) 


where imposed midplane stresses have not been included. As in Donnell's 
equation, equation (Cl) takes no account of inertia forces in the longi- 
tudinal and circumferential directions. 


Substitution of w and 1 for the unstiffened cylinder (eqs. (l), 
(3), and (If)) into equation (Cl) gives 


- a 2) + 

P R 

where the term 


Pi 

P 


c ± 

c 




has been 


effects of Sezawa (viscous) damping 


- F(M-k, v,n) - i -JL ek = 0 (C2) 

|v| 


added to include qualitatively the 
and where, now. 
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Equation (C 2 ) is precisely equation (8); hence, the flutter boundary in 
the supersonic range for the empty, unstressed cylinder is still defined 
by 



(C 4 ) 


but with n given by the more complicated function, equation (C 3 ) • As 
before, the critical boundary corresponds to the minimum value of M 
as v and n are varied. Setting 8 i) 2 /c*v = 0 yields 


-10 

v 


+ 3n 2 v® + 


2a 1 *- + 2(3 - 2|i)n 2 


12(1 - p2) 



v 6 


n c 


2rr 


2(5 - 2ji)n 2 + 4(2 - |x> - 


12(1 - H 2 ) 



3n^(n 2 - l) 2 v 2 - n^(n 2 - l) 2 = 0 (C5) 

Careful examination of equation (C 5 ) indicates that, for small values of 
h/R, there are two minimums of 2 2 with an intervening maximum. Computa- 
tions have been made for an aluminum cylinder at sea level by letting 
n = 2, 3, and 4 ; both the mininrums of fi 2 were checked and, through equa- 
tion (C 4 ) , the boundaries corresponding to the lower of these were deter- 
mined. These results are shown by the dashed curves in figure 4 ; in 
addition, the corresponding wave lengths are given in figure 5. 
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APPENDIX D 


IMPORTANCE OF DAMPING IN UNSTIFFENED CYLINDER 


It is of interest to examine some of the results of an analysis of 
the empty, unstiffened cylinder -which takes no account of damping (as 
opposed to the analysis in the body of this paper -which treats a cylinder 
in the limit as damping approaches zero) . The equilibrium condition, 
equation ( 9 ), reduces to 


Ls — _ ft 2 ) - F(M-k,v,n) = 0 (Dl) 

p R ' ' 


and its constituent functions vary as in figure 2 except that now 
Im(L) = 0 throughout the range M - 1 < k < M + 1. 

A typical Nyquist diagram (resulting when k traverses the path 
of fig. 3) is shown in the following sketch which corresponds to the 
particular conditions v > 0, 0<M-l<fl<M, Re C^)]s-M-l > 

and Re(SL/Sk)j 5 _j i _^ L < 0; 









42 


NACA TN .3638 


(The infinitesimal, coimterelockvri.se semicircles correspond to infinites in- 
itial, counterclockwise semicircles traversed by k to exclude zeros of L 
on the real axis from the lower half -plane.) These particular conditions 
are illustrated because they show clearly the possible existence of a 
stable condition with ReCL)^.^^ > 0, a circumstance found impossible 

with damping present. The reason for the difference, of course, is that, 
for k > 0, the damping term is a negative imaginary quantity and even 
the smallest amount of damping shifts all the portion^of the diagram 
corresponding to the range M-l<k<M+l below the imaginary 
axis . 


_ Note that, for the case of zero damping, satisfaction of equilibrium 
(L =0) with k real is not always sufficient for definition of the 
boundary; for the conditions in the cited example, the roots of H - 0 
do not pass directly from the upper into the lower half of the k-plane as 
the other parameters are varied, but linger awhile on the real axis. 
Instability occurs when a root leaves the real axlB and enters the lower 
half of the k-plane. For the cited example, this definition corresponds 
to the simultaneous solution of the equations 


L = 0 


8l 

_ sr 

8k 


0 


1 




) 


(D2) 


in the range M - 1 < k < M. 

Another significant difference in the two analyses is that, when 
damping is ignored, the instability in the subsonic region is low-frequency 
flutter and no purely static instability of the cylinder is possible. Sim- 
ilar differences would be found for the fluid-filled cylinder if the 
damping were ignored. The fact that the two assumptions, e — >0 and 
e ~ 0, produce such dis continuously different results testifies to the 
importance of the presence of damping in the unstiffened cylinder. 
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(a) 


Subsonic flow (M < l). 




(b) Supersonic flow (M > l). 

Figure 2.- lypical variations with k of functions comprising L. ..(See 

eq. (9).) 
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Figure 4 . - Stability boundaries for empty, infinitely long, unstiffened 
aluminum; cylinder at sea level with no applied membrane stress. 
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Figure 6.- Topical variation of F-j. with k. 
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Figure 7 *- Critical divergence "boundaries for infinitely long, unstressed, 
ring-stiffened aluminum cylinder at sea level. 




